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AFGL-TR -84-0195 

Inference  of  Equatorial  Field-Line- Integrated 
Electron  Density  Values  Dalag  Vhlatlara 

D.  Anderson  1 

Air  Foree  Geophysics  Lab.,  Hanscom  Air  Foroa  Base,  Mi  0173Jf 

P.  M.  Kintner  aad  M.  C.  Kelley 
Scbool  of  Elactrlcal  Engineering,  Cornell  Oniveralty,  Ithaca,  KT  1*653 

Abstract.  The  nighttime  electron  density  integrated  along  a  aagnetlc  field 
line  at  very  small  L- values  (L  -  1.06)  la  Inferred  by  coopering  whistler 
dispersions  aeasured  from  a  sounding  rocket  with  aodel  ionospheric 
calculations.  At  a  local  tine  of  0500  LT,  the  electron  density  in  the  F-layer 
valley  was  found  to  be  about  1  x  103/cm3.  Ve  suggest  that  this  technique  can 
be  applied  to  earlier  tines  in  the  local  evening  to  determine  ionospheric 
conditions  which  benefit  the  growth  of  low  latitude  please  Instabilities. 

Introduction 

During  the  post  sunset  period  the  equatorial  Ionosphere  is  known  to  be 
turbulent  over  scale  sizes  from  tens  of  kilometers  to  tens  of  centineters.  The 
turbulence  (equatorial  spread-F)  is  believed  to  originate  at  long  wavelengths 
through  an  Interchange  process,  such  as  the  layleigh-Taylor  instability  or  the 
EzB  drift  instability.  In  brief,  a  plasma  depletion  develops  at  the  bottem  of 
the  F-layer  which  rises  like  a  bubble  through  the  F-reglon.  Many  of  the  long 
wavelength  properties  can  be  simulated  using  two  dimensional  models  (Ossakow, 
1981)  which  neglect  variations  along  the  field  line,  nonetheless,  the  prooess 
Is  known  to  be  three  dimensional.  ALTAIF  radar  measurements  show  dramatically 
that  the  plana  depletions  fully  extend  along  magnetic  lines  (Tsunods,  I960). 
In  addition,  the  field  line  Integrated  paraaeters  such  as  Pedersen  oonduotlvlty 
and  electron  oontent  are  Important  in  determining  the  instability  growth  ratea 
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(Anderson  and  Haarandel,  1979).  On  the  otbar  band.  It  la  very  difficult  to  • 

■eaaure  tbe  electron  density  below  tbe  dense  F-reglon.  In  this  paper  ve  apply 
an  old  technique  (measuring  whistler  dispersion)  to  validate  nodel  calculations 
of  equatorial  field  line  Integrated  electron  densities.  Ve  also  suggest  that  • 

neasurlng  whistler  dispersion  at  appropriate  local  tines  nay  be  useful  in 
determining  whether  or  not  anblent  Ionospheric  conditions  exist  which  benefit 
the  growth  of  low  latitude  plaaaa  Instabilities.  • 

VLF  waves  below  the  electron  gyrofrequency  propagate  dlsperalvely.  For 
oo editions  at  tbe  equatorial  Ionosphere  the  local  dispersion  la  proportional  to 
tbe  square  root  of  electron  density  and  tbe  total  dispersion  Is  proportional  to  % 

the  Integrated  square  root  of  electron  density  along  the  group  path.  Our 
approach  will  be  to  compare  observed  whistler  dispersion  characteristics  to 
three  dlnenslonal  models  of  the  equatorial  Ionosphere  (Anderson,  1981)  »nd  to  ip 

show  that  ve  not  only  can  verify  the  nodel  but  oan  additionally  determine  the 
eleotron  density  within  the  F-layer  valley.  Unfortunately  the  local  tine  of 

«* 

these  particular  neasureaents  Is  at  0500  LT  which  limits  their  inportanoe  to  9 

spread-F  theories.  However  the  technique  does  not  sppear  to  be  sensitive  to 
local  tine  and  It  could  easily  be  applied  In  the  post  sunset  period. 

Experiment  Description 


A  Taurus- Tomahawk  Bounding  rocket  (3*. 010)  was  launched  from  Punta  Lobos,  _ 

Peru  on  March  21,  1983  et  9:53:30  UT  to  a  saxlmum  altitude  of  827  km.  Its  ^ 
primary  mission  was  to  atudy  tbe  critical  velocity  effect  using  a  barium  abaped  ,~] 
charge.  In  addition  to  tbe  cbenleal  release  the  payload  contained  a  TLF 

electric  field  experiment  which  sensed  potential  fluctuations  between  two - 

spheres  separated  by  3  esters.  Tbe  launch  point  was  about  20  km  south  of  the  Codes 


magnetic  dip  equator  end  the  payload  traveled  to  the  vest.  To  e  good]. 
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approximation  the  payload  was  elosa  to  tbs  dip  aquator  during  tba  ant  Ira 
flight. 

A  sketch  of  tba  ralatlon  batwean  tba  payload  and  tba  aquatorlal  magnetic 
flald  la  abown  In  Figure  1.  The  payload  la  approxlaataly  oantarad  along  tba 
flald  line  Intersecting  tbe  payload.  low  auppoaa  that  a  lightning  atroka 
occurs  at  one  end  of  tba  field  line  which  initiates  a  whistler.  If  wa  sake  tba 
important  assiaptlon  that  tbe  whistler  propagates  parallel  to  tbe  magnetic 
field,  then  wa  expect  tba  rocket  to  see  a  series  of  whistlers,  each  with 
progressively  more  dispersion.  Tbe  first  whistler  trace  will  have  dispersion 
associated  with  tbe  square  root  of  electron  density  integrated  along  1/2  of  tba 
field  line,  tbe  second  trace  will  be  associated  with  3/2  of  a  field  line,  tbe 
third  trace  will  be  associated  with  S/2  of  a  field  line,  etc.  Ve  should 
emphasize  that  there  is  no  a  priori  reason  to  expect  whistler  propagation 
parallel  to  tbe  sagnetlo  field  but  tbe  observation  of  a  multiple  bop  whistler 
showing  dispersion  with  tbe  ratios  1/2,  3/2,  5/2,  etc.  would  be  strong  evldenoe 
in  favor  of  parallel  or  ducted  propagation.  Tbe  propagation  of  low  latitude 
whistlers  has  been  reviewed  by  Hayakawa  and  Tanka  (1978)  who  concluded  that 
multiple-hep  whistler  trains  are  rarely  observed  free  ground  stations  for 
nighttime  low  latitude  conditions  although  they  say  be  sore  cossonly  seen  free 
in  altu  observations. 

An  example  of  one  of  nine  sultlple-bop  whistlers  observed  during  tbe  10 
minute  flight  la  shown  in  Figure  2.  This  example  oocurred  at  365  km  altitude 
during  tbe  upleg.  Tbe  whistler  clearly  shows  five  traoes  and  a  sixth  is 
observable  in  tbe  raw  date.  Tbe  dispersions  for  tbe  first  five  traces  are  8.6, 
16 ,  30,  37,  58  seel/?  which  nearly  Batches  tbe  expected  ratio  1/2,  3/2,  5/2, 
7/2,  9/2.  Tbe  dispersion  along  tbe  field  line  say  be  calculated  by  dividing 
tbe  expected  ratio  into  tbe  observed  dispersion  and  then  averaging.  Tbe  result 
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for  the  dispersion  along  ooe  full  length  of  the  field  11m  la  10.8  see^2. 
This  procedure  baa  been  applied  to  8  other  multiple  hop  whistlers  and  the 
estimated  dispersions  are  shown  in  figure  5  as  x'a.  At  425  km  altitude, 
corresponding  to  L  *  1.065,  the  typical  dispersion  is  11  secV2.  This  value 
agrees  quite  veil  vith  ground  station  observations  during  the  0000-0700  LT 
period  ( Hay aka v a  and  Tanaka,  1978). 

For  path  lengths  greater  than  1/2,  the  whistler  trace  in  Figure  2  has  a 
Hissing  band  between  325  Hz  and  425  Hz.  Ve  believe  that  the  slsslng  band  was 
caused  by  attenuation  near  the  two  Ion  crossover  frequency  Osltb  and  Brloe, 
1964)  for  a  mixture  of  B+  and  0*  or  heavier  ions.  The  two  ion  crossover 
frequency  is  given  by  fx  *  fCH«-(*W*  ♦  Hr+/16»  for  0*  where  ♦  %•>  r  H. 
Since  we  expect  that  Hg+  «  Nq+,  the  aeoond  tens  can  be  neglected  and  fx  s 
»0+/NfaH-f.  The  H+  cyclotron  frequency  at  425  ka  altitude  on  the  dip  equator  is 
340  Hz.  This  sets  the  low  frequency  boundary  of  the  Hissing  band  and  suggests 
that  the  local  B+  plasma  density  la  about  5f  of  the  total  density.  As  the 
whistler  propagates  downward  and  avay  free  the  dip  equator,  the  crossover 
frequency  Increases  until  the  whistler  reaches  the  foot  of  the  field  line.  The 
H+  cyclotron  frequency  at  the  foot  of  the  field  line  (125  ka  altitude)  is  451 
Hz  which  also  suggests  an  R*  contribution  to  the  total  density  of  about  5f . 


Model  Calculations 

The  aeasured  whistler  dispersion  yields  an  electron  density  averaged  along 
the  propagation  path.  To  Interpret  this  result  we  have  employed  an  ionospheric 
model.  To  oaleulate  electron  densities  as  a  function  of  altitude,  latitude  and 
local  time,  the  tlae-dependent  ion  (0+)  continuity  equation  ia  solved 
numerically.  This  equation  la  given  by 

♦  *  •  ("i^i)  •  -  h 
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vb«r«  Nj_  ( *  M«)  la  tha  Ion  density;  Pi,  tb«  Ion  production  rata;  Li,  tba  ion 
loaa  rata;  and  V*,  tba  Ion  tranaport  velocity.  Sol  vine  aquation  (1)  at  low 
latltudea  nacaaaltataa  transforming  the  Independent  coordiataa  r,  6  and  Mo  a 
coordiate  system  parallel  and  perpendicular  to  £  (Andereon,  1973a, b).  Tba  aat 
of  coefflcienta  for  the  ion  continuity  aquation  ia  obtained  from  aodels  of  tbe 
neutral  composition,  ion  and  electron  teaperaturea  and  production,  loaa  and 
diffusion  rates  as  well  aa  £  z  £  drift  and  neutral  wlnda.  Briefly,  tbe  aodels 
are  as  follows: 

1)  Tbe  MS  IS  neutral  ataospberlc  aodel  (Hedln  at  al.,  1977)  la  used  to 
calculate  H2,  O2,  and  0  densities  and  tbe  neutral  taaperature  as  a  function  of 
altitude,  latitude  and  local  tlae. 

2)  Production,  loss  and  diffusion  rates  are  alailar  to  those  used  by 
Anderson  (1973b).  For  the  pbotolonization  coefficient  at  tba  top  of  the 
ataospbere,  P  ,  a  value  of  6.0  z  10"? sec”1  la  cboaen  to  represent  March  20-21 , 
1983  conditions.  Electron  and  ion  teaperaturea  are  obtained  from  equlnoz, 
aoderate  solar  activity  measurements  at  Jlcaaarca  (March,  1969)* 

3)  Inclusion  of  the  vertical  £  z  £  drift  velocity  first  observed  by 
Woodman  (1970)  and  studied  in  detail  by  Fejer  at  al.  (1979)  la  essential  in 
producing  and  maintaining  electron  density  distributions  observed  near  tbe 
magnetic  equator.  Figure  3  displays  an  £  z  £  drift  aodel  appropriate  for  solar 
maximum  conditions.  For  tbe  March  20-21  simulation  tbe  vertical  drift  velocity 
pictured  in  Figure  3  waa  multiplied  by  0.8  between  1600-2000  LT  and  by  1*25 
between  0200-0500  LT. 

A)  Tbe  asavaed  meridional  ooaponent  of  tbe  neutral  wind,  taken  from 
Anderson  and  Clobucbar  (1983)  la  poleward  during  tbe  day,  reaches  a  maximum 
velocity  ef  70  a/eec  at  1700  LT  at  ±  85°  lat.  and  la  sy»etrlc  about  tbe 
geographic  equator.  The  vlad  becomes  equatorward  at  1915  LT  and  reaches  a 


■axiaua  velocity  of  85  a/sec  at  2300  LT.  Tbe  tonal  wind  ooaponent  is  neglected 
because  tbe  aagnetlc  field  declination  la  aull  at  Jleaaaroa  longitudes. 

Calculated  electron  density  values  at  a  function  of  altitude  at  tbe 
■agnatic  equator  are  Illustrated  In  Figure  4  and  eoapared  to  ground 
observations.  In  tbe  early  aornlng  (0200  LT)  tbe  Jlcaaarca  radar  observatory 
■ade  a  series  of  density  aeasureBents  as  a  function  of  altitude  which  ah  owed 
tbe  F-reglon  peak  at  about  350  ka  altitude  eoapared  to  tbe  aodel  prediction  of 
300  ka.  nonetheless  tbe  overall  profile  and  absolute  values  agree  well  with 
tbe  aodel.  At  0500  LT,  tbe  tine  of  tbe  sounding  rocket  flight,  tbe  only 
Ionospheric  eleotron  density  aeasuraaent  available  was  the  value  of  foFj  frm 
Buancayo.  Tbe  altitude  of  f^  was  270  ka  In  agreeaent  with  the  aodel  while 
tbe  aodel  electron  density  was  slightly  saaller  than  the  lonosonde 
aaasureaents.  On  board  tbe  payload  was  a  Langs ulr  probe  sensitive  to  tbe 
relative  electron  density;  It  Indicated  that  f0F2  w*  between  280  ka  and  260 
ka.  Tbe  electron  density  aodel  at  0500  LT  Is  used  in  the  calculations  that 
follow. 

Tbe  tine-dependent  continuity  equation  Is  solved  awerleally  along  a 
nuaber  of  geoaagnetlc  field  lines  In  order  to  ealculate  electron  densities  up 
to  an  altitude  of  1000  ka  over  tbe  aagnetlc  equator,  it  each  point  on  tbe 
field  line  where  electron  density  is  calculated  we  also  oelculate  tbe  looal 
plaaaa  frequency,  fp,  given  by  ^/60.5  *  i«(el/a^)j  f#,  tbe  electron 
gyrofrequency,  and  tbe  whistler  group  velocity  given  by 


integrated  along  the  field  line  fro*  125  k*  altitude  in  the  northern  hemisphere 
to  125  km  altitude  in  the  southern. 


or 


v/f  T  «  D  -  |  /  ds  assuming  that  f#  »  f . 
vfe 


D  la  called  Disperiaon  In  units  (see1/2). 

The  calculated  dispersion  la  coopered  with  the  observed  dispersion  In 
Figure  5.  Ve  assume  two  different  values  for  the  electron  density  In  the 
F-layer  valley  1  z  103  el/ cm3  and  2  x  10*  el/o*3  which  produce  two  different 
curves  for  the  dispersion  as  a  function  of  rocket  altitude.  The  assumption  of 
an  F-layer  valley  electron  density  Of  1  x  103  el/ cm3  produoes  values  which 
agree  well  with  observed  dispersions.  There  la  a  suggestion  from  the  data 
above  350  km  that  a  slightly  smaller  value  (20  %)  of  electron  density  In  the  F- 
layer  valley  would  be  appropriate. 

Conclusions 


Oslng  whistler  dispersion  aeaauraments  we  have  verified  an  Ionospheric 
model  of  electron  density  at  0500  LT  and  set  the  value  of  electron  density 
within  the  F-layer  valley  to  be  slightly  smaller  than  1  x  103/cm3  for  magnetic 
field  lines  which  Intersect  the  equator  at  higher  altitudes.  The  results  are 
consistent  with  other  In  situ  measurements  which  Imply  electron  densities  in 
the  F-layer  valley  at  the  magnetic  equator  of  the  order  103  at  earlier  local 
times.  (MoClure  et  al. ,  1977;  Morse  et  al.,  1977,  larelsl  and  S* us* ex evict, 
1981). 

it  earlier  looal  times  (1900-2300)  there  are  Indications  that  the  electron 
density  at  F-reglon  valley  altitudes  (180-J50  la)  man  he  as  high  as  2  -  3  * 
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10*el/cm3.  These  density  values  were  measured  in- situ  by  rocket-borne  disc 
probes  end  a  Bennett  Ion  ease  epectroaeter  during  the  Brazilian  Ionospheric 
Modification  Experiment  (BIKE)  conducted  at  Ratal,  Brazil  in  Septeaber,  19B2 
(Marclal  and  Baerendel,  1983)*  When  these  measured  valley  densities  mere 
Included  In  flux-tube  Integrated  Pedersen  conductivity  and  electron  content 
calculations,  It  was  found  that  there  was  a  significant  Increase  In  Instability 
growth  times,  suggesting  that  the  lack  of  observed,  naturally-occurring 
equatorial  plasmas  depletions  may  be  related  to  the  relatively  high  eleotron 
densities  within  the  valley  region  (Anderson,  1983)* 

To  determine  whether  or  not  vhlstler  observations  might  be  a  way  of 
•remotely"  sensing  the  electron  density  values  within  the  valley  region  during 
the  1900  to  2300  LT  period,  we  have  calculated  dispersion  values  as  a  function 
of  altitude  at  the  magnetic  equator  at  2000  LT  Instead  of  0500  LT.  The  same 
two  valley  electron  density  values  were  chosen,  1  x  103el/oa3  and  2  x 
10*01/0*3,  and  the  same  Input  parwetera  discussed  above  were  Incorporated. 
Figure  6  presents  the  calculated  dispersion  profiles  for  these  two  cases  as 
well  as  the  calculated  electron  density  profile  at  the  magnetic  equator.  The 
post-sunset  enhancement  In  upward  JA  drift  Is  responsible  for  lifting  the  F 
layer  so  that  the  valley  region  at  the  equator  Is  below  350  km  at  this  local 
time. 

At  the  bottom  of  the  F- layer  (350  km)  the  difference  in  dispersion  Is 
appreciable.  Assuming  a  valley  density  of  1  z  103el/om3  gives  rise  to  a 
dispersion  of  2  sec^/2  while  an  assused  value  of  2  x  10*el/em3  yields  a 
dispersion  of  6. A  eecl/2,  with  Increasing  altitude  the  dispersion  values 
increase  because  both  the  electron  density  and  the  length  of  the  geomagnetic 
field  line  are  Increasing.  The  difference  between  the  two  dispersion  values, 
however,  decreases  because  the  valley  densities  oontrlbute  leas  and  less  to  the 
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total  field  line- Integrated  density.  It  appears  therefore,  that  a  rocket 
launch  near  the  sagnetlc  equator  In  tba  post-sunset  tlae  period  aeasurlng  both 
YLF  wares  and  in- situ  electron  denaltlaa  would  be  able  to  determine  •remotely* 
▼alley  densities  associated  with  the  field  lines  it  Is  intersecting.  Such  a 
capability  would  help  establish  Instability  growth- rates  and  determine  whether 
In  fact  high  valley  densities  are  a  detrlaent  to  equatorial  bubble  formation. 
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Figure  Captions 

Figure  1.  A  schematic  diagram  showing  the  relation  between  the  sounding  rocket 
payload  and  a  whistler  propagating  parallel  to  the  magnetic  field.  The 
payload  would  observe  whistler  dispersions  associated  with  1/2,  3/2, 
5/2,... of  a  full  magnetic  field  line  length. 

Figure  2.  An  example  of  one  of  nine  multiple-hop  whistlers  observed  by  the 
sounding  rocket  payload  34.010. 

Figure  3.  Vertical  ionospheric  drift  velocities  at  the  magnetic  equator  as  m 
function  of  local  time. 

Figure  4.  The  model  Ionospheric  electron  densities  are  compared  to  J10 
measurements  at  0200  LT  and  to  Buancayo  tjfi  measurements  at  0500  LT. 

Figure  5.  The  calculated  dispersion  for  the  model  Ionosphere  using  two  assisted 
values  for  the  F  layer  valley  (2  x  10*/o*3  and  1  x  103/em3)  compared  to 
the  dispersions  measured  on  the  sounding  rocket  payload  34.010. 

Figure  6.  Calculated  whistler  dispersion  profiles  for  equatorial  Ionospheric 
conditions  typical  of  2000  LT  assuming  two  different  values  for  the 
F-layer  valley  electron  density  (lx  103/cm3  and  2  x  1oVcm3). 
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